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Abstract. We study possible Lorentz violations by means of gamma-ray bursts 
(GRB) with special focus on the Large Array Telescope (LAT) of GLAST, where 
we concentrate on models with linear corrections to the speed of light. We simulate 
bursts with gtobssim and introduce a Lorentz violating term in the arrival times of 
the photons. We further perturb these arrival times and energies with a Gaussian 
distribution corresponding to the time resp. energy resolution of GLAST. We 
then vary the photon flux in gtobssim in order to derive a relation between the 
photon number and the standard deviation of the Lorentz violating term. We 
conclude with the fact that our maximum likelihood method as first developed in 
[l] is able to make a statement whether Nature breaks the Lorentz symmetry if 
the number of bursts with known redshifts is of the order of 100. However, the 
systematic errors caused by unknown mechanisms for photon emission are not 
considered here despite the fact that these errors should be the main obstacle to 
detecting Lorentz violations. 
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1. Introduction 

The space-time structure as described by general relativity is a classic one in the sense 
that it is smooth even at very small distances and very large energies. However, there is 
an agreement that new phenomena caused by quantum gravity (QG) should affect this 
smoothness at a distance of the order of the Planck length Ip — \J hG / ~ 1.6 x 10"^'^ 
cm or equivalently the Planck mass Mp = y^hc/G « 1.2 x 10^^ GeV/c^. Despite 
the fact that we still lack a full quantum theory of gravity we can make some 
predictions derived from phenomenological quantum gravity. One of these predictions 
is a distortion of the photon dispersion relation [51 [31 13] 



where E denotes the photon energy, p its momentum, a a model-dependent 
dimensionless parameter of order unity and c the speed of light in vacuo. This 
distortion leads to an energy-dependent velocity of light given by 



A more general relation with an explicit Lorentz symmetry breaking could be used 
[9], however there are stringent experimental limits on first order in Ip modifications 
to dispersion relations in a lorentz breaking scenario, see e.g. [llj or the constraint 
from the Crab [TUl [H]. On the other hand Lorentz invariance could be seen as a 
consequence of the foamy structure of space-time caused by quantum fluctuations on 
short times and distance scales fl3l [HI [151 [16] . Another approach, called deformed 
(or doubly) special relativity (DSR) [17l [181 [El EHl |21] , introduces a second constant 
scale (the Planck length Ip or equivalently the Planck energy Ep) that is observer 
independent. This yields a non-linear Casimir operator of DSR which may result 
in an energy-dependent speed of light, depending on the chosen model. However, a 
common feature of all DSR models is that the speed of light does not depend on the 
helicity. 

One reason why present experiments have not been able until now to either rule 
out or confirm a violation of the Lorentz symmetry is the fact that the correction to the 
velocity of light in Eq. ^ is suppressed by the Planck energy. However, despite this 
small effect it was pointed out that one powerful way to look for an energy-dependent 
velocity of light is given by gamma-ray bursts (GRB) [4l[5l[6l[7l[8l[22l|23l|24l|25l|26]. 
GRBs are the most luminous electromagnetic events occurring in the universe since 
the Big Bang. They can last from a few milliseconds to minutes, with a typical 
duration of a couple of seconds [571 HH]- However, the main problem we face when 
looking for QG effects in GRB signals is our ignorance of the internal physical processes 
which are at the origin of the photon emission. Photons with different energies may 
emanate from different mechanisms, thus further complicating the comparison between 
events. Despite these difficulties several studies were able to put strong constraints 
on a Lorentz invariance violation with GRBs ,221 [301 [33 [32] ■ In a recent work |33] 
a preferred range for the linear QG mass scale of Mqg ~ 0.4 x 10^^ GeV could even 
be found by studying a flare observed by MAGIC. This result has a sensitivity that 
probes, for the first time, the Planck mass range. 

In this work we shall mainly consider the first-order correction to the speed of 
light with the conservative bound a — 1. In other words we assume that a Lorentz 



E^ = p^c^ + a— + 0{E^/El), 



(1) 




(2) 
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symmetry breaking starts to become important only at the Planck scale Mqq = Mp, 
resulting in the following relation for the speed of light: 

The dependence of the speed of light on the energy as described by Eq. ^ causes 
a difference in the arrival times between two photons emitted at the same time, but 
with different energies. Assuming a flat universe described by the ACDM model (for 
different models see |34]), the integrated time delay between two photons with an 
energy difference /S.E is given by [35] 

A^ = ±i/-i^^ r^^±i^=dz', (4) 
° Mpc^ Jo Vf^A + + z'f 

where flm = 0.27, JIa = 0.73 and Hq = 71 km Mpc~^. At the end of the paper we 
will briefly discuss the case of a quadratic correction where the integrated time delay 
is given by 

The paper is organized as follows. In section[2]we describe the relevant properties 
of GLAST and in section 13.11 and section 13.21 our method to create a photon list with 
arrival times and energies. Section is devoted to explaining our results, section !?^ 
to giving an estimate of GLAST to quadratic corrections to the speed of light and 
section 15751 to studying non-FRED distributions. We then conclude with section [H 



2. GLAST 



The Gamma Ray Large Area Space Telescope (GLAST) is a space-based gamma-ray 
telescope designed to explore the high-energy universe. It includes two instruments: 
the Large Array Telescope (LAT), which is an imaging gamma-ray detector which 
detects photons with energy from about 30 MeV to 300 GeV, and the Gamma-ray 
Burst Monitor (GBM) that consists of 14 scintillation detectors which detect photons 
with an energy between 8 keV and 30 MeV. The LAT has a very large field of view that 
allows it to see about 20% of the sky at any time. Despite the fact that it will cover 
the entire sky every three hours we shall assume that the instrument response changes 
on timescales longer than a typical burst duration. On the other hand, GLAST can 
be pointed as needed when a bright GRB is detected by either LAT or GBM so that 
it will detect around 200 GRBs each year. The energy resolution ranges from 20% at 
30 MeV to about 7% at 1 GeV, as can be seen in figure [T] [36]. The time resolution of 
an event should be around 10 /zs with a dead time shorter than 100/is. In summary, 
LAT will have superior area, angular resolution, field of view, time resolution and 
deadtime. This will at least provide an advance of a factor 30 in sensitivity compared 
to previous missions. 

The GLAST Science Support Center (GSSC) provides analysis tools freely 
available to the scientific community on their homepage [37] . As we are interested in 
simulating the detection of GRBs by GLAST we mainly used the tool called gtobssim 
which is a software that generates photon events from astrophysical sources with the 
instrument response functions of GLAST. Because we will only study the measurement 
of possible Lorentz violation we assumed that LAT pointed in the same direction as 
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the burst. Further information on the effects of QG on LAT GRBs can be found in 

e.g. El mi- 




Figure 1. Energy resolution as a function of the energy for the LAT |36) . The 
energy uncertainty at 30 MeV is about 17% before going down to a couple of 
percent in the GeV range. Most of the photons LAT detects have an energy 
around ~ 100 MeV with an uncertainty of about 10%. 



3. Creation of a photon list 

3.1. Analysis method 

We shah use the same method we used to study INTEGRAL GRBs by modehng bursts 
with a Fast Raise and Exponential Decay (FRED) distribution [1]. We parameterize 
this distribution with four parameters: / — f{ti, Ef, P, R, D, k, h) as shown in figure[21 
We suppose that a photon i came from the probabihty density function / at time ti 
with energy Ei . We use the method of maximum hkehhood which consists in finding 
the set of values P, R, D, k and h that maximizes the joint probability distribution 
of all data, given by 

T{P,R,D,K,h)=\[f{t,,E,-P,R,D,K,h) (6) 

i 

together with the constraint 
ti 

dt' f{t',E,-P,R,D,n,h) = l, (7) 

D 

where J- is the likelihood function and the integral runs between to and ti as shown in 
figure [21 The specifics of this distribution does not play a significant role in our study 
as we are only interested in seeking possible violations of the Lorentz symmetry. In 
order to further improve our model we also simulated an isotropic background with a 
spectrum following an exponential decay with exponent 2.1. However, this background 
does not play a significant role because of the following two reasons. The first one is 
that LAT shall only be able to detect photons with energies above ^ 30 MeV, where 
such events are rare. The second one is the fact that we are only interested in GRBs, 
e.g. events of short duration. 
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Figure 2. Sketch of a typical liglit curve of a GRB for a given energy interval. 
The curve is parameterized by four parameters: R the duration of the rise, h the 
height above the background, D the decay time for exp(— t/D) and k describes the 
magnitude of the dependence on the energy of the distribution /, t 

max — -\- K - ^ 

where P is the time when the intensity reaches a maximum and E is the photon 
energy. The dashed line shows a distribution for another energy interval that is 
shifted by an amount of At = k ■ Ai? sketching the shift in time due to quantum 
gravitational effects. 




Figure 3. Left panel: Example of a simulated LAT GRB with gtobssim with a 
total photon number of 74. The burst starts at to = 5 s and has a decay time of 
the order of 10 seconds. Right panel: Spectrum of the same GRB as simulated by 
gtobssim for LAT. The first detected photons have an energy of roughly 30 MeV. 
We choose the high energy exponent (3 = 2 (see Eq. ||8)) and a background with 
exponent 7 = 2.1 



3.2. Spectrum of the GRBs 



As described in the previous section we also need to simulate the energy of the photons. 
Normally a typical energy distribution of GRBs follows the pattern of the so-called 
Band function [S^ given by the following equation: 



Ne{E)=A 



A 



E 



100 keV 



(a - P)Eq 



cxp 



E 

To 



(a - l3)Eo > E, 



100 keV 



E 



100 keV 



exp(/? — a), 
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(« - mo < E, (8) 

where a is the low-energy exponent, (3 the high-energy one and Eq the break energy. 
As LAT starts measuring at 30 MeV and the break energy is around 500 keV we shall 
only be interested in the high-energy behavior of the Band function. However, LAT 
will open a new window on the spectrum where little is known, therefore there is no 
certainty whether the Band function is still valid throughout the energy range of LAT. 

In order to make the simulations as realistic as possible we used gtobssim. As 
described in the previous section we simulated a GRB with photons following a FRED 
distribution together with an isotropic background. Since we are only interested in 
the detection of a possible Lorentz violation we used the same FRED distribution for 
all GRBs, i.e. with a raise time of the order of a second and a decay time around 10 
seconds (see figure [3]) . In order to get a feeling of the uncertainty we varied the flux 
of the burst and only selected the bursts which got about the same number of hits in 
the detector. We then introduced a Lorentz violating term with the following relation 
(see figure [5]): 

At = KxE-y, (9) 

where is the energy of the detected photon in MeV and At is the time delay in 
seconds caused by quantum gravity given by the relation ti = + Ai. The parameter 
K describes the effect of a Lorentz violation and is given by the following relation: 



QG<: 



1 + z' 



/^3 



-.dz' =: 



M, 



(10) 



where Mqg is the mass scale where the Lorentz symmerty breaks down. Henceforth 
we shall take the most conservative scale and set Mqg = Mp « 1.2 x 10^^ GeV. 

As can be seen in figure |4] typical values for k. lie between 10~^ and 10~^ s/MeV. 
In the Monte Carlo simulations we used a value of k = 4 x 10~^ corresponding to a 
redshift oi z — 1. 




Figure 4. The Lorentz violation parameter function of the redshift. As 

we shall see below, a value of k = 10~* s/MeV could even be measured without 
requiring a bright burst (see table [TJ. 



We studied two models: one simple one with only the time delays caused by 
quantum gravity and a more realistic one where we perturbed the energy and the 
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arrival time according to the scheme shown in figure [5] in order to take into account 
the energy and time resolution of the LAT. For a given event at time with energy 
Eq obtained with gtobssim we read the energy uncertainty I from figure [T] and 
perturbed this energy with a Gaussian with maximum at Eq and standard deviation 
■ The next step is to take the corrected time ti and perturb it with a Gaussian with 
maximum at t = and standard deviation at = 10 fxs so that we get the perturbed 
arrival time t2- 
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Figure 5. Setup of the simulation process. The software gtobssim gives as 
output the arrival time to and the energy Eq of a photon. The arrival time is 
then corrected by means of the relation l[9) from to to ti (box QG). The simpler 
case with only QG effects is shown by the arrow number 1. In the second case, the 
energy of the photon was perturbed from Eg to Ei with a Gaussian distribution 
with standard deviation according to figure [l] (box ER). Furthermore, the arrival 
time was perturbed a second time from ti to t2 with a Gaussian distribution with 
standard deviation of 10 /is corresponding to the time resolution of LAT (box 
TR). The photon with arrival time t2 and energy Ei is then given to the fitter 
(box FIT) to get a value for the QG corresponding to k. 



3.3. Results 

In the previous section we explained how we constructed a burst with photons and 
perturbed their energies and arrival times to account for Lorentz violation and finite 
energy and time resolution of LAT. The question is now whether it is possible to get 
back the value of k parameterizing the Lorentz violating term despite the perturbation 
of both the energy and arrival time. In [1 we found an exponential dependence 
between the standard deviation ct^ and the number of detected photons N with an 
exponent of -0.617. However, we only perturbed the arrival time and not the energy. 

The final step is to search for R which minimizes the likelihood For a single 
burst we scanned through a large range of k in order to find the global (and not just 
a local) minimum of the likelihood. We performed 100 Monte Carlo simulations for 
different luminosities of the bust, i.e. different numbers of total events N (see table[T]). 

Figure [6] shows a comparison between the standard deviation of the perturbed 
and unperturbed system. The fit of the unperturbed system is given by 

a, = (2.5 X 10-2) . Ar-0.93 (11) 

and the fit of the perturbed system by 

cr« = (6.0 X 10-2) . (12) 

The results described above allows us to simulate a full set of GRBs with different 
redshifts z and time shifts k. Figure [7] shows a set of 100 GRBs with redshifts 
between and 10, where we used the results from [ID] and [321 i'^ order to get a 
rough estimate of the expected number of GRBs detected by GLAST as a function of 
the number of detected photons and the redshift. Each value for k has been perturbed 



GLAST and Lorentz violation 



8 



unperturbed (1.) perturbed (2.) 



"NT 

JN 


fl r„— 1 —21 
nux [s m J 


K [s/MeVJ 


CTk [s/MeVJ 


K [s/MeVJ 


(7^ [s/MeVJ 


20 


10 


1.2 X lO"'' 


2.0 X 10-3 


3.9 X 10-4 


1.7 X 10-3 


50 


28 


9.5 X lO"'"^ 


9.5 X 10-^ 


6.4 X 10-5 


7.8 X 10-4 


75 
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100 


60 
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4.3 X 10-4 
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3.7 X 10-4 
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96 


4.2 X IQ-s 


2.2 X 10-4 


4.6 X 10-5 


3.3 X 10-4 
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128 


5.3 X 10-6 


2.1 X 10-4 


3.7 X 10-5 


2.4 X 10-4 


500 


330 


2.9 X 10-5 


6.6 X 10-5 


1.4 X 10-5 


6.6 X 10-5 


1000 


660 


2.3 X IQ-'^ 


5.1 X 10-5 


2.9 X 10-5 


4.3 X 10-5 



Table 1. Results of 100 Monte Carlo simulations. The first column shows the 
number of photons detected by LAT, the second one the approximate flux in 
m~^s~'^. The third column shows the mean value of R for the unperturbed system 
(1.) and the fourth one its standard deviation. The last two columns are the same 
as the third and fourth ones, except for the fact that the arrival times and energies 
were perturbed with a Gaussian distribution (see section [3.21 1. 



10"' 




10 



10 10 



Figure 6. Comparison between the standard deviations (t„ of the unperturbed 
and perturbed system. The vertical crosses show the values for the unperturbed 
system together with the solid line as its fit. The inclined crosses show the values 
for the system with perturbed arrival times and energies. 



with a Gaussian with a standard deviation given by Eq. (jl2p . Considering a hnear 
approximation to quantum gravitational effects we have the relation |30[ 131] : 

K = al{z) + h{l + z), (13) 

where a and b arc fitted coefficients. The constant b parameterizes time lags in the 
frame of the source caused by unknown internal processes of the GRBs, a describes 
the expected Lorentz violating effects through 

TT-l 

ath = —2-^ « 2.56 X 10-5 s/MeV. (14) 

A fit of the 100 GRBs shown in figure [7] gives 

n = (2.6 ± 3.1) X 10-5 . 1{z) - (1.9 ± 2.5) x 10-5 • (1 + z). (15) 
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Figure 7. Simulation of 100 GRBs wliere k has been perturbed with a Gaussian 
according to the relation I I12I I. The error bars show the Icr deviation given by 
Eq. II12I I. The distribution of the GRBs and the luminosity has been computed 
using results from I40|[39l where a histogram of the number of GRBs as a function 
of either the redshift or the luminosity is given. The solid line shows the fit given 
by Eq. JTsJl. Note that no systematic errors have been incorporated. 



We see that the Lorentz violating term is well compatible with the theoretical value 
ath while the second term is less than one a off the input values 6th = (see table [5]). 
Considering the fact that LAT should see around 500 GRBs each year it might be 
tantalizing to conclude that only after a couple of months the question whether the 
Lorentz symmetry is broken might be answered. However, we would like to stress the 
fact that we did not consider systematic errors caused by the lack of knowledge of the 
internal processes leading to photon emissions. As these systematic errors were the 
main problem past works [H [30l [3TJ [32] had to deal with, this oversimplified analysis 
must be taken with precaution. Moreover, a known redshift of the bursts is also 
needed, thus reducing the number of usable bursts. 

N a [s/MeV] b [s/MeV] 



10 (1.7 ± 5.6) X 10-5 (-1.8 ± 3.3) x 10"^ 

20 (2.6 ± 3.6) X 10-5 (-1.8 ± 2.8) x 10"^ 

50 (2.6 ± 3.2) X 10-5 (-1.9 ± 2.6) x 10-^ 

100 (2.6 ± 3.1) X 10-5 (-1.9 ± 2.5) x 10-5 



Table 2. Results of the fit 11311 for a variable number of bursts. The first 
column shows the number of bursts considered for the fit, the second one shows 
the value of a describing the Lorentz violating effects and the third one the value 
of b parameterizing the time lags caused by unknown internal processes. The 
theoretical value of a is given by Eq. The error on a and b only reduces 

considerably between N = 10 and = 20. 



3.4- Quadradic corrections 

In this subsection we are concerned with quadratic corrections given by Eq. ([5]). Such 
corrections may be interesting in view of the fact that previous works have already 
put stringent constraints on linear corrections to the speed of light [lOl [TTJ [32l [33] . To 
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get a rough estimate of the sensitivity of GLAST to quadratic corrections we must 
first get a bound on the time difference that should be detectable by GLAST. With 
Eq. (HD we see that for cr^ - 5 x 10"^ s/MeV (see table [J) and l\E ^ 10^ MeV we 
get a bound on the time difference of (XAt 5 x 10~^ s. Solving Eq. ([5]) for the mass 
scale and inserting these results we get a mass scale of 

Mquadratic > 2 X 10^ GcV. (16) 

for z = 1 up to which GLAST should be sensitive for a single burst. This bound may 
probably be raised by a couple of orders of magnitude with better statistics. However, 
we do not think that GLAST will be able to detect a quadratic correction if the QG 
scale is at the Planck scale. 

The above estimate may seem very crude. We therefore checked it with results 
obtained in the literature and found a good agreement between this estimate and a 
more rigorous statistical analysis. For example, taking Ml > 7 x 10^^ GeV obtained 
in Eq. (36) in jJOj together with /\E - 100 keV and cta* 5 x 10"^ g foj. bATSE we 
get a bound on the quadratic correction of Mq > 0.8 x 10^ GeV, which corresponds 
more or less to the result Mq > 3 x 10^ GeV obtained in [3^ with a rigorous analysis. 

3.5. Other distributions 

Until now we only considered FRED distributions for GRBs. However, despite the 
fact that these distributions may describe a large number of GRBs, we also have to 
consider other burst shapes in order to get an estimate for the validity of our likelihood 
method with non-FRED distributions. We studied two other shapes, a linearly raising 
and falling distribution and a double-peak distribution (see figure IS]). We then applied 
the likelihood method to these peaks with 200 photons and compared the results with 
table [1] for = 200. For the linear peak (see left panel in figure [8|) we found a 
value K = 1.7 X 10"^ s/MeV and a standard deviation of ctk = 1-6 x 10^^ s/MeV. 
Comparing this value with the perturbed system for N — 200 (see table [T]) we see 
that the likelihood method gives a slightly better results for the linear peak. We also 
studied the double peak (see fight panel in figure[8]) and found a value oiR = 2.0 x 10~^ 
s/MeV with a standard deviation of = 1-6 x lO^"'. Thus the results of both cases 
lead to the conclusion that our likelihood method is also able to deal with non-FRED 
distributions. 

We fixed the duration time of the simulated GRBs to about 20 seconds, which 
seems quite arbitrary. However, we expect the likelihood to get better as the duration 
of the burst gets shorter (as long as the total number of photons remains constant). 
The reason is that, while k is independent of the duration, the ratios k/R and k/D 
(see figure[2]) are inverse proportional to ti—t^. We checked this claim for bursts with 
a duration time of about 6 seconds, 200 photons and FRED-distributed. We found a 
value K = 2.6 X 10~^ s/MeV with a standard deviation of CTk = 3.7 x 10~^ s/MeV. 
Comparing this value with table [T] we see that the likelihood is able to recover the 
Lorentz violating term with a precision of about one order of magnitude better. On 
the other hand, the convergence of the likelihood decreases with the duration time of 
the bursts, thus narrowing the study of short bursts to the brighter ones. 

4. Conclusion 

As can be seen from table [T] values of the Lorentz violation parameter k of the order 
of ~ 10-5 s/MeV could theoretically be measurable. In [39l gQl |4T] the BATSE 
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Figure 8. Left panel: Example of a simulated LAT GRB with a linear raise and 
linear decay. Right panel: Example of a simulated LAT GRB with a double-peak 
shape. 



catalog was used to build up the statistics for LAT. The authors made the assumption 
of extrapolating the BATSE observations to LAT energies in order to get a rough 
estimate on the number of detected bursts each year as a function of the number of 
photons per burst and the energy threshold (see Table 1 in [iQ]). If this assumption is 
correct we may get a couple of bursts each year that could answer the question whether 
the Lorentz symmetry is violated. Nevertheless, as explained above, the systematics 
should considerably wash out the signal from a Lorentz violating term so that a lot 
of bursts with known redshifts should be needed to get a tight bound on a Lorentz 
violating term. Even without systematics a number of bursts of the order of 100 still 
yields an uncertainty of the order of the expected Planck corrections, as can be seen 
from Eq. (fTS]) . 

In this work we had to make a couple of simplifications. For example we had 
to assume that the BAND function is still valid at LAT energies. In other words, 
GLAST will be able to measure the high-energy index at energies out of reach of past 
and present experiments. Furthermore, we studied the effect of a linear correction to 
the speed of light and assumed a simultaneous emission of the photons of different 
energies from the source. This last assumption is the most problematic one because in 
[T] we obtained results for k with even different signs for the same burst (see also e.g. 
[30l[3Tl|42] where the same kind of problems arises). Therefore, the unknown emission 
mechanisms for photon emission should be the major problem in seeking Lorentz 
violation. This is the reason why either a lot of bursts or a better understanding of 
the burst mechanism is needed before we can settle down the question whether the 
Lorentz symmetry is violated in Nature. 
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